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Abstract —In this paper, a mathematical framework for the 
analysis of average symbol error prohahility (ASEP) in millimeter 
wave (mmWave) cellular networks with Poisson Point Process 
(PPP) distrihuted base stations (BSs) is developed using tools from 
stochastic geometry. The distinguishing features of mm Wave com¬ 
munications such as directional beamforming and having different 
path loss laws for line-of-slght (LOS) and non-llne-of-sight (NLOS) 
links are incorporated in the average error probability analysis. 
First, average pairwise error probability (APEP) expression is 
obtained by averaging pairwise error probability (PEP) over 
fading and random shortest distance from mobile user (MU) to 
its serving BS. Subsequently, average symbol error probability 
is approximated from APEP nsing the nearest neighbor (NN) 
approximation. ASEP is analyzed for different antenna gains and 
base station densities. Finally, the effect of beamforming alignment 
errors on ASEP is investigated to get insight on more realistic cases. 

1. Introduction 

Demand for mobile data and use of smart phones have been 
increasing very rapidly in recent years. According to the UMTS 
traffic forecasts, 1000 fold increase in mobile data traffic is 
predicted by the year 2020 HI. In another estimate, more than 
50 billion devices may be connected wirelessly by 2020 which 
may cause a capacity crisis m. To meet this increasing demand, 
moving to new frequency bands becomes a necessity. Therefore, 
the large available bandwidth at millimeter wave (mmWave) 
frequency bands, between 30 and 300 GHz, becomes a good 
candidate for the fifth generation (5G) cellular networks and 
has attracted considerable attention recently IJl, ID, Q, Q. 

Evaluating the system performance of mmWave cellular 
networks is a crucial task in order to understand the network 
behavior. There are several recent studies which analyze the 
coverage probability and average rate in mmWave cellular net¬ 
works using results from stochastic geometry and the theory of 
point processes for different base station(BS)-user associations 
Q, la, HI. Stochastic geometry is a commonly used powerful 
mathematical tool to evaluate the average network performance 
of spatially distributed nodes ifTOI . Poisson point process (PPP) 
is a widely used model in wireless networks in general and 
in cellular networks in particular due to its analytical tractabil- 
ity. However, average error probability in PPP-based cellular 
networks has only been barely analyzed in the literature. For 
instance, there is work focusing on the computation of average 
symbol error probability (ASEP) in the presence of Poisson field 
interferers (see e.g., HD). However, none of them are applicable 
to cellular networks since the BS to mobile user (MU) cell 
association is generally not considered. In ifT^ . a mathematical 
framework to compute the ASEP in cellular networks, where the 


BS locations are modeled as independent homogeneous PPPs, 
is established for the first time. Their approach is based on 
the shortest BS-to-MU distance cell association criterion, which 
guarantees that the interfering BSs are located farther than the 
serving BS, so it is applicable to cellular networks. However, 
to the best of our knowledge, average error probability analysis 
has not been conducted for mmWave cellular networks yet. 

In this paper, we follow a similar approach as in na to 
develop a mathematical framework for the computation of 
ASEP in downlink mmWave cellular networks. First, average 
pairwise error probability (APEP) is calculated by averaging 
PEP over fading and random shortest distance from MU to 
serving BS. Then, ASEP is found using the nearest neighbor 
(NN) approximation. The main contribution of this paper is the 
combination of the characteristic features of mmWave commu¬ 
nications with the proposed mathematical framework in HD- 
One distinguishing feature of mmWave cellular communication 
is the directional beamforming at the transmitter and receiver, 
which provides an array gain to mitigate the effect of path loss. 
In this paper, sectored directional antenna model is used to find 
the effective antenna gain similar to Q, il, HI- First, perfect 
beam alignment is assumed between the MU and the serving 
BS. Then, the effect of beamsteering errors is investigated. 
Another distinct feature of mmWave communication is that a 
BS can be in the line-of-sight (LOS) or in the non-line-of-sight 
(NLOS) of MU and different path loss laws are applied for LOS 
and NLOS links. Here, we adopt the equivalent LOS ball model 
in Cl to determine whether a BS is LOS or NLOS. 

H. System Model 

In this section, we introduce our system model for the down¬ 
link mmWave cellular network consisting of BSs distributed 
according to some homogeneous PPP 4/ of density A in the 
Euclidean plane. Without loss of generality, we consider that 
a typical MU is located at the origin. A shortest distance cell 
criterion is assumed, i.e., MU is served by the nearest BS which 
is denoted by BSq. The distance from the ith BS to the MU 
is denoted by for i G 4/. Thus, the distance between the 
MU and serving BS (BSq) is rg which is a random variable 
(RV) with PDF froiO = 27rA^exp{— lfl5l . The set of 
interfering BSs i G 4*—BSg is still a homogeneous PPP, denoted 
by vl/(\o)^ according to the Slivnyak-Mecke’s Theorem ifTOl . 
We assume that all the interfering BSs are transmitting in the 
same frequency band as the serving BS (full frequency reuse), 
therefore has density A as well. 


We have the following two assumptions in the construction 
of the system model. 

Assumption 1: Antenna arrays at both the BSs and MU 
are used to perform directional beamforming such that the 
main lobe is directed towards the dominant propagation path 
while smaller sidelobes direct energy in other directions. For 
tractability in the analysis, antenna arrays are approximated by 
a sectored antenna model, in which the array gains are assumed 
to be equal to a constant M for all angles in the main lobe and 
another smaller constant m in the side lobe Qa. The MU and 
serving BS, BSq, are assumed to have perfect beam alignment 
and therefore have an antenna gain of MM. Also, the beam 
direction of the MU and each interfering BS can be modeled as 
a uniform random variable on [0,27r]. Therefore, the effective 
antenna gain is a discrete RV described by 

( MM with prob. pmm = (^)^ 

Gi = < Mm with prob.pMm = 2^2^ ^ ( 1 ) 

[ mm with pmh. p^m = 

where 6 is the beam width of the main lobe. 

Assumption 2: A BS can be either LOS or NLOS BS to the 
MU according to the LOS probability function p(r) which 
is the probability that a link of length r is LOS. Using field 
measurements and stochastic blockage models, p{r) can be 
formulated as where decay rate /3 depends on the building 
parameter and density El. LOS probability function p{r) 
can be approximated by a step function in order to simplify 
the analysis. In this approach, the irregular geometry of the 
LOS region is replaced with its equivalent LOS ball model with 
radius Rb Q. A BS is a LOS BS to the MU if it is inside the 
ball, otherwise it is a NLOS BS. Different path loss laws are 
applied to LOS and NLOS links. Thus, the path-loss exponent 
on each interfering link can be expressed as follows: 


I OL if r < i?B 
ajq if r > Rb 


( 2 ) 


where ar and un are the LOS and NLOS path-loss exponents, 
respectively. 

By combining these two assumptions with the described 
network model above, the received signal at the MU can be 
written as. 


y = y GoEor^ “^/iqSo + y/O^r-hiSi+n. (3) 

' r: ' 

--,-' 

Lgg 


where x is the signal arriving from the serving BS to MU, Jagg 
is the aggregate network interference, and n is the Gaussian 
distributed noise component with zero mean and variance Nq. 
Moreover, Go is the effective antenna gain of the BSq-MU link 
and it is assumed to be equal to MM, Eq is the BSs’ transmit- 
energy per transmission, aL is the LOS path-loss exponent 
of the BSq-MU link, sq = aoexp{j0o} is the information 
symbol transmitted by BSq with amplitude oq and phase 9q, 
ho = |/io| exp {j(/)o} is the fading coefficient in the BSq- 
MU link where |ft.oP is an exponential RV with parameter 
CTo = Ej/iop = 1 and the phase (fo is a uniformly distributed 
RV in the range [0, 27r). A similar notation is used for /agg, but 
note that the effective antenna gain Gi and path loss exponent 


ai are different for different interfering links as described in ([T]i 
and dU, respectively. Also, considering phase modulation, we 
assume that oq = = 1 for i € 

At the MU, an interference-unaware maximum-likelihood 
(ML) demodulator is used as in El, which can be formulated 
as O 

So = argmin{i:)(so) = |y - \/GoEor;f°^^hoSol^j- (4) 

So 

Inserting Q into (|4]l and neglecting some irrelevant constants 
after algebraic manipulations, we can express the decision 
metric as 

D(so)cxry“^GoL;o|A,,,n/io|2 

-f 2r(7“^yG^Re{(/agg -f n)|/io|exp{-j(/)o}A*^^}, 

(5) 

where = sq - sq, /agg = VGiM)r~°‘* hiSi and 

n ^ CAf{0, No). Since the effective antenna gain between the 
MU and each interfering BS is modeled as an independent 
RV, we can employ the thinning property of PPP to split the 
aggregate network interference /agg into 6 independent PPPs as 
follows Iflrl : 

T _ / jMM _|_ jMM \ j_ { jMm _|_ jMm \ \ ( jmm , jmm \ 

V'^^LOS ■^^NLOS/ V'^^LOS '^^NLOS/ v'^^LOS “^^NLOS/ 

= X! (-^^LOS + ^^NLOs)> (6) 

GG{MM,Mrn,rnm} 

where each interfering BS is either a LOS or NLOS BS and the 
superscripts represent the discrete random antenna gain defined 
in O. According to the thinning theorem, each independent 
PPP has a density of Xpa where pc is given in ([T]i for each 
antenna gain G £ {MM, Mm, mm}. 

III. Average Error Probability Analysis 

In this section, we investigate the error performance of a 
downlink mmWave cellular network. The first step in obtaining 
an approximation of the average error probability is to com¬ 
pute the pairwise error probability (PEP) associated with the 
transmitted symbols. Hence, initially we derive an expression 
for PEP, conditioned on fading gain (|/io|) and random shortest 
distance of the MU-serving BS link (rp), in terms of the char¬ 
acteristic function (CP) of the aggregate network interference 
and the noise. A closed-form expression is determined for the 
CP of the aggregate network interference for PPP distributed 
BSs. Then, APEP is computed by averaging the conditional PEP 
over fading and the position of the serving BS. Pinally, ASEP 
is approximated from APEP using the NN approximation. 

A. Derivation of Pairwise Error Probability 

The PEP is defined as the probability of deciding in favor of 
So when actually sq is transmitted. It is assumed that these two 
symbols are the only two symbols in the signal-constellation, 
and therefore decision is made strictly between these two 
symbols. Using the decision metric in (|5]l, PEP conditioned on 
|/io| and ro can be computed as 

P{so so||/io|,To} = P{D{so = So) < D{so = sq)}. (7) 

When So = sq, Ag ^ = sq — sq becomes zero by definition. As 
a result, D(so = sq) is zero. Let U = /agg + n. Note that U is 
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a circularly symmetric RV. Thus, U exp {j0o}^g{^* s} = U 
US). Thus, PEP can be computed as follows: 

P{so So||^o|,fo} = P{D{so = So) < 0} 

= p!^Re{U\ho\exp{jMK,s} < 

= p|Re{t/}<-^^^|A,,,||/io|} 

= Fu.^[-^^^\AsM\ho\^ (8) 

where Fu^ denotes the CDF of the RV C/rs = Re{t/}. 

Gil-Pelaez inversion theorem can be employed to compute 
the CDF Fu^ by using the CF of Urs, $[4 (w), as follows 

ES): 


further express it using similar steps as in im as follows: 

(w) 


= exp{XpcTr{R% - rl)[-l + E^^^ri{exp{jw\/GEQr^ “^-zj}]} 

= exp |2ApG7rE2,{ J {exp{jw^/GEQr~°“^Zi} - l)}ridri| 
= exp |2ApG7r J ($z(\/ GEowr~°“^) - l)ri(iri| 

= exp |2ApG7’‘ J ($o(\/GF;ok|r““^) - l)ridri| 


exp < —2XPCTT 


(yG^|zu|)2/“^ 1 _ 


CtL 


/GEo\w\R^ 


F/o‘L + 1 


dt 


Fus^iu) 



^l_ir 

2 TV Jo w 

Im{(cos(wu)—jisin(wu))(Re{'l?GRe(u’)})+iIm{^t/Re(l^)})}^^ 

w 


sin(wu)$u(w) , 
—^i— -aw. 
w 


(9) 


where the last equation follows from the fact that the CF 
<i)GRc('tt') is ^ real function, i.e., Im{$GRe('tf)} = 0 ^nd 
Re{<i)Gp,^(r(;)} = $GRe(ti'), and C7(.) is a circularly symmetric 
RV, i.e., ^Urc{w) = ^u{w)- 


(13) 

where (a) follows from having a PDF of 2ri/(i?^ — rg) 
if ro < Ti < Rb and zero otherwise, (b) originates from the 
definition of the CF, (c) follows from the fact that z has a 
spherically symmetric (SS) PDF and its CF is also SS, i.e., 
‘hz(re) = ^’o('tf) for some $o(-)’ (d) is obtained by applying a 
change of variables with t = \/GEo\w\r~‘^^. 

$o(f) can be found using the properties of an SS distribution: 


By inserting (|9ll into ([ 8 ]l, PEP can be written as = Ea;,{cos(fa:i)} + j Ej^,{sin(fyj)} 


P{so^So\\ho\,ro}=J-- f sin w 

2 Wo V 2rg J 

In ( fTOb . PEP is computed using the CF of the RV U. Since U 
is the summation of two independent RVs, Jagg and n, ^uiw) 
is equal to the product of the CFs of these two RVs, i.e., 
T>g(w) = $ 4 gg(w)^'n(w)- ^'n(w) = exp{-W(7Vo/4)} is the 
CF of a Gaussian RV ll20l and is calculated in the next 

subsection. 


B. Characteristic Function of the Aggregate Interference Rgg 

Since /agg is the sum of six independent PPPs as seen in (| 6 l), 
using stochastic geometry, its CF can be written as 

^/.gg(r^)= n (w), (11) 

^LOS ^NLOS 

where $jg (w) and <i)jG {w) are the CFs of FOS and NFOS 

’I'LOS ^los 

components with antenna gain G. 


0 

= Ea,,{cos(fa:i)} (14) 

where Xi = Re{zi}, yi = lm{zi} and the second term is zero 
because sin is an odd-symmetric function. 

By inserting the result in (fT4l l into (fOl l and taking the 
expectation operator Ea:;!.} outside, the integral inside the 
exponential function can be calculated as shown at the top of 
next page in ( fTSl l using the result from ll22l Eq. (3.771.4)], 
where pFg is the generalized hypergeometric function. Then, by 
inserting the result of the integral in (fTsT i into (fTJt and applying 
similar steps as in |[T2), we obtain the closed-form expression 
for $jG (w) in (fT^ . 

^LOS 

Similarly, a closed-form expression for the CF of the interfer¬ 
ence due to NFOS BSs, 4 >jg {w), can be obtained by changing 

^LOS 

the boundaries of the integral and replacing with a at in 
(fTST l. More specifically, since NFOS BSs lie outside of the ball, 
integral is evaluated from Rb to infinity and the expression for 
‘hjG (w) is found as shown in (fTTl i on the next page. 

^LOS 

Finally, a closed-form expression for the CF of the aggregate 
network interference, can be obtained by inserting 

equations ( fTSl l and ( fTTl l into (fTTT i. 


Fet Zi = Sihi = \hi\exp{j{0i +fi)}. The interference 
due to a FOS component with a generic antenna gain G can 
be written as = Ejg^los VG'£'or““^Zi. Then, its CF 

<i)jG {w) = E{exp{jiy/^^^^}} can be obtained using the same 
steps as those in ED and can be expressed as in (fT^ given on 
the next page. By using the Taylor series expansion for the 
exponential function, one can rewrite the equation in (fl^ and 


C. Average Pairwise Error Probability 

In this section, APEP is computed by averaging PEP. Aver¬ 
aging can be performed by taking the integral of the conditional 
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T,; = 


= i: - rm^PO’^(K - ’■g)] (E..„{exp{i^'/G^r.-«».}))'-'. 

k^O 

1 - cos{tx^) i _ cos(te,) 


( 12 ) 


^ (v^kl) 


l2laL + l 
— ‘Ijalj 


0 
2 

S ' 0 


l2laL + l 


7?2 -7.2 

XId / , 


' "V «l’2’ cl’ 4r2“- V ' 'V aL’2’^ a,.’ 


(15) 


(w) = exp{XpGTT{rl - Rl)} 


I 2 7^ / 1 1 1 1 1 Gi?o|ti'p0'0 \ 2 /I 111 1 Gi?o|wPcro 

X exp<{ -ApG7rro2^2 ( —;- 7 ^ 177 ;— ) + ApG7ri?|2^2 ( - 7 :,-;:7,1-;-- 


2’ ai’2’ ai’ Arl^^ 


2’ aL’2’‘ az,’ 


$zG (w) = exp{ApG7’‘-R|} exp < -XpgttR% 2 F 2 


_1_ l_.l ^ _1 ^ GEo\w\'^ao 

2’ a7v’2’ ttiv’ 4i?|“" 


(16) 

(17) 


PEP over |ft,o| and tq as follows: 

Eavg{so -)■ So} = E|^j,|_ro{E{so so||/io|,ro}} 

= E|/*ol,'-o ^ (^^^^^\As,s\\ho\w'^w~^A>u{w)dw^ 

= \-\ E-o {e|.oI {sin } ^c/(«;)| 


Also, let denote the number of nearest neighbors of sq 

having distance Amin- Now, ASEP can be approximated as 

ASEP « kd^i^P^yg{Amin}- (23) 

where = 2 when modulation order (7) is greater than 2, 
w~^dw, and A^m = 2sin(7r/7) assuming multilevel phase shift keying 
(18) (MPSK) modulation. 


where (a) follows from the fact that $g(w) depends only on 
ro not |/io|. Hence, the expectation over |/io| can be computed 
in closed-form by employing the PDE of Rayleigh distribution 
and calculating the resulting integral as m 


El 


l^ol 


sin 


sin 


y/GoEo 

( VGqEq 


\ 2r“- 


|As_s||/io|?« 

As,s||/lo|w 


CTo (To 


-^/GqEqGq 

4r“- 


|As,s|wexp<^ - 


GqEqGq I 

Te^ 


|A,.,|W . (19) 


By substituting ( fT^ and the PDF of ro (i.e., froiO) (flST l. 
APEP can be expressed as follows: 


Eavgjso —>■ So} ^ 2 ~ 




V GoEoao 
45“r 


As.s|wJexp < — 


GoEgao . 


27rA^exp{— ttA^ }4>(7(w)ui d^dw. 


|2 2 
W 


( 20 ) 


Finally, substituting ^ij{w) into ( l20b and after some algebraic 
manipulations, APEP can be rewritten as shown at the top of 
next page in dSTT i. where we dehne SNR = i?o(To/4 and 


2F2ix,y) 


(I 11 1 GSNRiwpA 

V2’"y’2’^"y’ )■ 


( 22 ) 


D. Average Symbol Error Probability 

In this section, we approximate ASEP from APEP in (ISTT l 
by using NN approximation. The advantage of the NN approx¬ 
imation is that it only depends on the minimum distance in 
the constellation and the number of nearest neighbors 12^ . In 
Section II, |As s| is dehned as the distance between the constel¬ 
lation points So and sq. Hence, we dehne A^m as the distance 
of So to its nearest neighbors, i.e., Amin = minso^^o lA^o.^ol- 


IV. ASEP IN THE Presence oe Beamsteering Errors 

In Section mil MU and the serving BS are assumed to 
be aligned perfectly and ASEP is calculated in the absence 
of beamsteering errors. However, in practice, it may not be 
easy to have perfect alignment. Therefore, in this section, we 
investigate the effect of beamforming alignment errors on ASEP. 
We employ an error model similar to that in ll24ll . Let |e| be the 
random absolute beamsteering error of the MU’s beam toward 
the serving BS with zero-mean and bounded absolute error 
Emax A tt. It is appropriate to consider the absolute beamsteering 
error due to symmetry in the gain G. The PDF of the effective 
antenna gain G with alignment error can be explicitly written 
as ID 

/G(g) = F|,| (^0 5{g-MM) 

+ ‘^F\A (Jj)S{g-Mm) 

+ ^(g-mTO), (24) 

where 6{-) is the Kronecker’s delta function, F\^\{x) is the 
CDF of the misalignment error and (l24l) follows from the 
dehnition of CDF, i.e., E'|e|(a;) = P{|e| < a:}. Assume that 
the error is distributed according to a Gaussian distribution, 
so absolute error |e| follows a half normal distribution and 
E’|£|(a:) = erf(a:/(-\/2(TBE)), where erf(-) denotes the error 
function. 

From (EB, it is clear that PEP depends on the effective 
antenna gain between the MU and the serving BS, and so does 
the ASEP. Thus, PEP can be calculated by averaging over the 
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r. < ~ 1 VGoSNR,^ i-oo .oo . SNR,^ ,2 o'! r r 

fkvg{so So} = - - y/^x- — -|As,s| J J exp|-^|2^^|A^,spw2| exp{-7rA^2}gxp{_y;2jVo/4} 

X [exp{XpcTT^'^}exp{-XpcTT^'^2F2{^,aL) + XpGnR%[2F2iRB,aL) - 2F2{RB,aN)]}f'^d^dw. (21) 

G^{MM ,Mm,mm} 




Fig. 1: ASEP as a function of the SNR = Eoaoji for different BS densities 
A (tti = 2.1, ajv = 4, M = 10 dB, m = -10 dB, 9 = 15, BPSK ) 


distribution of G, fois)’ follows: 

-Pavg{so —^ So} = E(5{Pavg{so —^ Sq; G}} 

^^avglso -)> So;g}/G(g)rfg 

= (F|,|(0/2))2Pavg{so^so;MM} 

+ 2(F|,|(0/2))F|,|(0/2)Pavg{so ^ so; Mm} 

+ .F|e|(6'/2)^Pavg{so-)■ so;mm}, (25) 

where we define F\^\{6/2) = 1 — F|g|(0/2). 

V. Numerical Results 

In this section, we provide numerical results to evaluate the 
error performance of a downlink mm Wave cellular network. In 
all figures, LOS and NLOS path loss exponents are ul = 2.1 
and aN = 4, respectively. In the non-mmW case, all BSs are 
assumed to be LOS and the path loss component is equal to 2.1. 
Also, the radius of the LOS ball Rb is assumed to be equal to 
141 meters similarly as in Q. 

First, we compare the performance of the mmWave network 
with that of the non-mmWave network (antennas are omnidi¬ 
rectional, and all BSs are LOS). In Fig. 1, ASEP versus SNR 
is plotted for different BS densities with BPSK modulation. 
As shown in Fig. 1, we have better error performance in the 
mm Wave scenario than in the non-mmWave one. Also, with the 
increasing BS density, ASEP is decreasing. 

Next, we plot the ASEP with different antenna main lobe 
gains and different BS densities. The numerical results in Pig. 
2 show that with increasing main lobe gain M, ASEP decreases 
significantly. Also, note that different combinations of main lobe 
gain and BS density, e.g. (M = 20dB, X = 10“^) and (M = 
lOdB, X = 10“"^) lead to the same error performance. Hence, 
the same error performance can be achieved by either decreasing 
BS density while increasing the main lobe gain, or vice versa. 

In Pig. 3, we also compare ASEP for different modulation 
orders assuming MPSK modulation. As the modulation order 



Pig. 2: ASEP as a function of the SNR = Eocro/4 for different antenna main 
lobe gains M and BS densities A (a^ = 2.1, ojv = 4, m = —10 dB, 9 = 15, 
BPSK ) 



0 


10 20 30 40 50 

SNR (dB) 


Pig. 3: ASEP as a function of the SNR = EqooH for different modulation 
orders Mo = 2.1, = 4, M = 20, dB m = -10 dB, 9 = 15, 

A = 10--*) 

increases, the minimum distance between the nearest neighbors 
decreases. Thus, as expected, error performance of the network 
gets worse with the increase in modulation size. 

Pinally, the effect of beamsteering errors on the error per¬ 
formance is analyzed in Fig. 4. ASEP versus SNR is plotted 
for different standard deviations of the alignment error. As 
can be seen, ASEP is getting worse with the degradation in 
the alignment angle, ctbe = 0 corresponds to the case with 
no alignment error and it has the best error performance as 
expected. Since ctbe = 2 has the same ASEP with ctbe = 0, we 
can infer that the alignment error until ctbe = 2 can be tolerated 
and ASEP increases significantly for ctbe > 2. 

VI. Conclusion 

In this paper, we have analyzed the average error perfor¬ 
mance of downlink mmWave cellular networks, incorporating 
the distinguishing features of mmWave communication into 
the average error probability analysis. Sectored antenna and 
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Fig. 4; ASEP as a function of the SNR = Boo’d/ 4 in the presence 
of beamsteering error for different standard deviations of alignment error 
{ol = 2.1, = 4, M = 20, dB m = -10 dB, = 15, A = 10“®, 

BPSK ) 

simplified ball-LOS models have been considered to simplify 
the analysis. Numerical results show that employing directional 
antennas improves the error performance. Also, we show that 
better ASEP values can be obtained by increasing BS density 
and main lobe gain. Investigating the effects of using different 
LOS probability functions instead of the simplified ball-LOS 
model, and incorporating more general fading models remain 
as future work. 
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